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ABSTRACT 
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been  added.  The  variation  of  the  transition  point  distance  with  the  shore  distance 
as  a  variable  was  observed.  In  additionjthe  isopressure  patterns  were  verified  by 
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I.  INTRODUCTION 


Sound  propagation  in  an  ocean  with  a  sloping  bottom  has  received  considerable* 
attention  during  recent  years.  Both  tHeoretlcIa^s  and  exp^imentalists  have  been 
involved  with  this  subject  [Ref  1—19].  They  have  created  a  number  of  acoustic 
models  to  predict  the  sound  field  within  a  wedge  shaped  fluid  overlying  a  penetrable 
or  a  rigid  fluid  bottom. 

The  acoustic  field  in  a  wedge  shaped  shallow  water  duct  with  ideal  boundary 
conditions  has  been  studied  by  Bradley  and  Hudimac[Ref.  1).  They  have  analyzed 
the  case  of  an  isospeed  duct  with  one  pressure  release  and  one  rigid  surface.  The 
theoretical  analysis  has  been  carried  out  in  both  image  theory  and  normal  mode 
theory.  The  relationship  between  the  two  has  been  explored  for  the  case  of  a  point 
source.  The  image  sum  was  transformed  via  the  Poisson  transformation  formula  to 
the  normal  mode  sum.  The  average  acoustic  field  w’as  then  compared  for  the  two 
approaches  with  good  results.  Also  the  two  theoretical  approaches  were  compared 
with  experimental  data  from  the  west  coast  of  Florida  with  fair  results. 

Kuznetsov  [Ref  2]  used  the  method  of  images  to  investigate  the  high-frequency 
field  in  a  wedge.  He  showed  that  the  method  of  images  is  better  suited  to  the  near 
field  where,  in  addition  to  the  normal  modes  of  the  discrete  spectrum,  the 
continuous  spectrum  plays  a  significant  role.  At  distances  from  the  source  greater 
than  ten  times  the  layer  thickness,  the  normal  mode  method  works  better.  These 
results  are  based  on  the  idealized  model  of  a  w’ave  channel  as  a  layer  of  water  having 
plane  boundaries  over  a  fluid  half  space.  The  sound  field  in  a  wedge  was  described 
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as  the  sum  of  the  image-source  fields  taking  account  of  the  coefficient  of  all 
reflections.  A  point  source  was  used  in  a  tiredge  with  perfectly  reflecting  boundaries. 

Kuznetsov  in  another  study  [Ref  3],  has  described  the  displacement  of  the 
normal  modes  which  takes  place  in  the  vertical  plane  and  causes  each  normal  mode 
to  emerge  from  the  wedge  into  the  half  space  at  a  certain  distance,  specific  to  the 
normal  mode,  from  the  vertex  line  of  the  wedge. 

In  1975  Graves,  Nagl,  and  Uberall  [Ref  6)  used  the  adiabatic  range  variation 
method,  derived  by  Pierce  and  Milder  [Ref  4,5],  to  perform  an  approximate 
separation  into  normal  modes  for  the  problem  of  an  under -ocean  channel  with 
gradual  range  dependence  of  medium  and  boundaries.  Their  technique  was 
illustrated  by  application  to  the  isovelocity  ocean  wedge  with  rigid  ocean  floor,  and 
compared  with  the  exact  solution  of  Bradley  and  Hudimac  [Ref  1).  A  good 
agreement  was  obtained  for  moderately  large  wedge  angles,  even  when  mode 
coupling  was  neglected. 

In  1978,  Coppens,  Sanders,  loannou  and  Kawamura  [Ref  7).  using  computer 
models,  predicted  the  acoustic  pressure  amplitude  and  phase  in  the  upslope  direction 
at  the  bottom  of  a  wedge  shaped  fluid  layer  overlaying  a  fast  fluid.  The  slopes  of 
the  wedges  studied  were  about  2.4®  or  2.7®.  One  of  the  models  had  infinite  source 
distance  w’hile  the  other  had  the  source  at  finite  distance.  The  program  outputs 
were  compared  with  the  experimental  results  using  as  variables  the  wave  number 
and  density.  There  were  significant  differences  between  the  simple  model  and  the 
experimental  results  because  of  the  experimental  limitations. 

The  results  of  shallow  water  acoustic  experiments  performed  off  the  coast  of 
Corpus  Cristi,  Texas,  have  been  compared  with  theoretical  computations  by  Rubano 
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(1980)  [Ref  8].  A  four-layer,  fluid,  normal— mode  model  was  used  to  predict  th^ 
group  velocity  dispersion  and  spatial  amplitude  distributions  of  the  first  and  second 
modes.  The  results  show  good  agreement  between  the  measured  and  theoretical 
values. 

In  the  same  year  (1980),  comparison  between  experiment  and  theory  was  made 
by  Jensen  and  Kuperman  [Ref  9].  They  used  the  parabolic  equation  model  to  study 
the  modal  cutoff  during  upslope  propagation  in  a  wedge-shaped  ocean.  They  found 
that  the  cutoff  is  not  abrupt  but  takes  place  over  a  finite  distance  which  essentially 
provides  an  aperture  for  radiation  of  a  beam  into  the  bottom.  Also  the  Gaussian 
representation  of  a  point  source  in  the  PE  model  has  given  very  good  results  in  the 
presence  of  a  bottom. 

The  modal  cutoff  during  upslope  propagation  has  also  been  studied  by  Pierce 
(1982)  [Ref  10].  Since  the  adiabatic  mode  theory  for  upslope  propagation  in  shallow 
water  breaks  down  when  the  depth  decreases  to  a  critical  value,  he  has  used 
matched  asymptotic  expansion  techniques  to  yield  the  acoustic  pressure  in  the 
transition  case.  The  derived  characteristic  critical-depth  transition  function  give 
results  which  are  in  accord  with  computations  reported  by  Jensen  and  Kuperman 
[Ref  9). 

The  same  subject,  the  modal  cutoff  during  upslope  propagation,  was  studied  by 
Arnold  and  Felsen  using  the  image  theory  as  a  convenient  starting  point  for 
collective  treatment  of  ray  fields  and  their  conversion  into  local  modes  [Ref  11]. 
The  numerical  comparisons  with  the  results  of  Jensen  and  Kuperman  [Ref  9]  showed 
good  agreement. 

The  sound  field  in  an  absorbing  fluid  substrate  underlying  a  wedge  shaped  fluid 
with  higher  sound  speed  has  been  studied  by  Coppens,  Sanders,  and  Humphries  in 


3 


1984  [Ref  12].  A  Green's  function  extension  of  the  image  theory  give  the 
opportunity  of  doing  all  the  calculations  on  a  computer.  Two  different 
approximation  methods  were  used.  The  end-point  approximation  gave  solutions 
valid  at  great  distances,  and  for  absorptions  representative  of  sedimentary 
materials.  The  saddle-point  approximation  has  given  good  results  under  the  apex, 
with  no  limit  in  the  absorption. 

In  1984  Baek  [Ref  13]  predicted  the  pressure  amplitude  and  phase,  everywhere 
within  the  wedge,  in  upslope  direction,  having  a  fast  bottom.  Three  cases  were 
studied,  pressure  release  bottom,  rigid  bottom,  and  penetrable  bottom. 

In  the  same  year  (1984)  LeSesne  [Ref  14]  studied  and  compared  two  computer 
programs  using  the  method  of  images  to  determine  the  pressure  and  phase 
distribution.  It  w'as  concluded  that  an  infinite  source  is  not  approached  until  the 
source  exceeds  six  hundred  dump  distances  from  the  apex.  All  measurements  were 
taken  directly  upslope  with  source  distance  equal  to  forty  dump  distances  and  the 
receiver  on  axis  with  the  source.  Comparison  between  the  two  programs  show 
similar  results. 

One  year  later,  in  1985  Borchardt  [Ref  15]  experimentally  measured  the  sound 
pressure  field  everywhere  within  isospeed  water  overlying  a  sloping  absorbing 
bottom.  The  experimental  data  were  compared  with  the  program  developed  by 
LeSesne  [Ref  14].  The  results  showed  that  the  program  successfully  predicts  the 
on— axis  upslope  sound  profile  in  a  wedge-shaped  medium.  Particularly, 
measurements  showed  good  agreement  for  R2  <  10  ,  /?  =  9.5°,  ^2/^2  ~ 

Y0<50  .  All  the  symbol  definitions  are  given  in  the  List  of  Symbols  table. 

A  shallow  water  wedge  with  a  slope  up  to  9°  has  been  experimentally  studied 
by  Tindle.  Hobaek,  and  Muir  in  1986  [Ref  16].  They  presented  experimental  results 
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for  downslope  propagation  and  they  defined  the  wedge  modes  as  propagating  normal 
modes  with  wavefronts  which  are  curved  into  arcs  of  circles  centered  on  the  wedge 
apex. 

In  1987  Kaswandi  [Ref  17],  studied  the  slow  bottom  case.  The  outputs  were 
taken  for  the  on-axis  case  with  wedge  angles  of  6®,  10®,  and  15®.  Three  kinds  of 
curves  were  observed,  depending  on  the  source  distance. 

At  about  the  same  time,  a  program  with  3-D  capabilities  was  used  by  Li  [Ref 
18]  to  predict  the  acoustic  pressure  distribution  along  the  interface  between  a 
wedge-shaped  fluid  layer  overlying  fast  or  slow  bottom.  It  was  concluded  that  the 
major  characteristics  of  the  output  graphs  were  insensitive  to  minor  variations  of 
the  bottom. 

In  1988  Doolittle  and  Tolstoy  presented  experimental  data  obtained  in  East 
Australian  Continental  Slope  consistent  with  the  theory  of  energetic  horizontal 
refraction  due  to  multiple  reflections  from  a  sloping  bottom  [Ref  19].  When  viewed 
from  above,  the  path  of  the  ray  appears  to  curve.  As  a  result,the  exact  position  of 
the  source  can  be  obtained  from  the  details  of  this  curvature. 

The  purpose  of  this  research  is  to  take  all  the  outputs  provided  by  the  theses  of 
Kaswandi  and  Li  and  try  to  build  a  physical  understanding  of  the  acoustic  field  from 
them.  Also  an  attempt  was  made  to  obtain  detailed  comparisons  between  the 
DSLOW  and  3LS4  programs  for  the  same  inputs,and  to  further  investigate  the 
transition  point. 
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II.  THEORY 


The  theoretical  background  of  this  r^earch  has  been  presented  by  Coppens. 
Sanders,  loannnou  and  Kawamura  [Ref  8].  What  follows,  in  this  chapter,  is  a 
summary  of  the  most  important  points  of  the  image  theory  and  its  application  to 
the  wedge-shaped  ocean  with  penetrable  bottom. 

A.  METHOD  OF  IMAGES 

The  assumptions  that  (1)  the  speed  of  sound  is  the  same  everywhere  within  the 
wedge  and  in  the  penetrable  bottom  and  (2)  the  upper  surface  is  a  pressure  release 
surface  makes  the  method  of  images  an  appropriate  approach  for  understanding  the 
sound  field  in  a  wedge-shaped  fluid.  The  sound  paths  are  straight  lines  and  the 
plane  wave  Rayleigh  reflection  coefficients,  discussed  by  Kinsler.  Frey.  Coppens 
and  Sanders  [Ref  14],  is  a  good  approximation  if  the  sound  source  is  not  too  close  to 
the  bottom. 

The  geometry  of  the  model  used  to  predict  the  pressure  amplitude  and  phase 
everywhere  within  the  wedge  is  shown  in  Figure  1.  For  the  wedge-shaped  duct, 
cylindrical  coordinates  are  used  with  the  shore  line  as  the  axis  (Figure  2).  The 
images  lie  on  a  circle  centered  on  the  shore  line. 

In  Figures  3  and  4,  the  relationship  between  the  receiver  and  the  nth  image  is 
shown.  The  angle  of  the  nth  image,  measured  from  the  bottom,  is  given  by: 

Bn  =  B'(n-1)+G,  for  n  odd 
6ln  =  B*n-G  for  n  even 
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According  to  Sommerfeld's  concept  of  an  extended  Riemann  surface,  the 
extension  of  the  $  coordinate  makes  the  number  of  images  infinite.  The  total 
pressure  and  phase  at  any  field  point  within  the  wedge  is  found  by  coherein 
summation  of  all  spherical  waves  radiated  from  the  infinite  number  of  images, 
the  number  of  reflections  associated  with  a  given  image  increases,  the  pressure 
contribution  of  this  image  decreases.  The  infinite  summation  is  approximated  be¬ 
taking  the  sum  over  N  images  where 

n  =  INTEGER  (ISO  /  P) 

Wedge  angles  of  the  form  r/n  have  been  used,  where  n  is  an  integer,  because 
diffraction  terms  were  eliminated  [Ref  2]. 

The  distance  Rn  between  the  receiver  and  the  nth  image  is  given  by: 


R 


n 


Ri2+  Ra^  -  2‘RiR2-cos(^„-D)  +  Yo* 


for  the  upper  group  of  images,  and 


Rn 


Ri2+  Ro^  -  2'R,R2-cos(<?„+D)  +  Yo- 


for  the  lower  group  images. 

The  angle  of  incidence  (|)nni  of  the  n***  image  on  the  m^**  plane  (  described  by  the 
angle  m  •  /? )  is 


7 


Sin(|)nm  — 


[Ri-sin(^n—  m,3)+  R2-sin(mi  -  D)] 

R^ 


The  reflection  coefficients  R(^n)  for  the  upper  family  of  images  and  R(^n')  for 
the  lower  family  of  images  for  a  plane  w’ave  are 


R(^nni) 


P\C\ 

~ 

'ML  r 

PiC2 


^nm 


^niB 


and 


R(^n'ni) 


PlCl 

'M2  ~ 

£l£i  + 
P2C2 


'I'n'.n 

^n'm 


where 


^nm  s 


y  1  -  (Ci/C2)Ws^(gnm) 
sin(^nii.) 


and 

q.  •  y  1  -  (Ci/C2)^COS^(gn'iB) 

^  n  m  = 

sin(^n'iD) 

The  resulting  complex  pressure  from  the  upper  family  of  images  is 
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N 

Pu=  ^ 
n=l 


1 


R(^n!r:) 

m=0 


and  the  resulting  complex  pressure  from  the  lower  family  of  images  is 


Pi  = 


N 


U=1 


M 

n  RiOn'J 

m=0 


The  sound  field  radiated  by  each  image  has  been  multiplied  by  the  plane  wave 
reflection  coefficients  corresponding  to  reflections  encountered  by  the  wave  as  it  goes 
from  the  source  to  the  field  point. 

The  total  complex  pressure  distribution  P  is  the  sum  of  the  Pu  and  Pi: 

P  =  Pu  +  Pi 


B.  NORMALIZATION 

All  distances  are  normalized  in  terms  of  the  characteristic  distance  X..  Tiii.'  is 
the  distance  measured  from  the  apex  along  the  wedge  interface  to  the  point  at  which 
the  lowest  mode  attains  cut  off  as  shown  in  figure  5  [Ref  20]. 

The  characteristic  distance  is  defined  as  the  distance  from  the  shore  at  which 
the  depth  H  is  given  by 


k-H 


It 
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For  a  fast  bottom,  the  critical  angle  6c  is 

=  COS^c 
C2 

since 

t 

kH  = - 

2sin^c 

tan/?  =  - 

X 

jr 

k  X  - - 

2  sinf^c]  tao;^ 

The  above  formula  gives  the  derivation  of  the  characteristic  distance  X  in  the 
case  of  a  fast  bottom  (Cj  <  Cq).  When  slow  bottom  case  is  studied  (Cj  >  c.,).  the 
characteristic  distance  can  not  be  defined  as  above.  A  useful  characteristic  distance 
for  slow  bottom  can  be  defined  by  inverting  Ci/C2 

k  X  - - ^ - 

2  tan[arcos(c2/Cj)]  tan/? 

The  advantage  of  working  with  normalized  distances  is  that  the  analysis  is 
frequency  independent. 

A  convenient  normalization  for  the  pressure  amplitude  is  accomplished  by 
multiplying  the  pressure  amplitude  at  one  meter  from  the  source  by  the  source 
distance  (Rl).  This  simplifies  the  interpretation  considerably,  especially  in  the  case 
of  large  source-apex  distance  compared  to  X. 
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In  the  course  of  this  research,  a  computer  program  has  been  created,  the  DEME 
program.  This  computer  program  provides  an  additional  normalization.  It 
normalizes  all  the  pressure  amplitudes  to  the  first  pressure  measurement  taken  down 
from  the  surface.  This  slope  normalization  makes  the  shapes  of  the  curves 
comparable,  even  if  they  have  been  obtained  at  different  distances. 
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III.  COMPUTER  PROGRAMS 


Two  programs  have  been  used  for  the  purpose  of  this  research.  The\-  are 
improved  and  more  documented  versions  of  the  programs  discussed  by  Kaswandi 
[Ref  17]  and  Li  [Ref  18].  They  use  the  same  model  but  have  different  output 
displays. 

A.  PROGRAM  DEME 

In  the  beginning  of  the  study  reported  in  this  thesis,  the  DSLOW  program  was 
obtained  [Ref  17].  DSLOW  is  itself  an  improved  version  of  the  WEDGE  and 
XSLOPE  programs  developed  by  Kawamura  and  loannou  [Ref  7],  and  LeSesne  [Ref 
14]  respectively. 

DSLOW  was  developed  to  run  on  the  IBM  3033  main  frame  of  the  Naval 
Postgraduate  School.  It  uses  the  method  of  images  to  predict  the  pressure  amplitude 
and  phase  within  the  wedge  fluid  overlying  a  slow  bottom. 

DEME,  created  from  DSLOW,  uses  the  same  model  to  predict  the  pressure 
amplitude,  and  phase  within  the  wedge  fluid  overlying  a  slow  or  a  bottom. 

The  output  is  given  in  a  table  format  or  in  a  graph  as  shown  in  the  Appendices. 
The  DISSPLA  graphics  package  and  the  EASYPLOT  program,  both  of  them 
available  on  the  IBM  3033.  have  been  used  for  the  graphical  outputs. 

During  this  research,  other  versions  of  the  original  DEME  program  have  been 
created  with  different  outputs.  All  of  them  have  been  stored  under  Professor 
Coppens'  account  in  the  IBM  3033  in  the  Naval  Postgraduate  School. 
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The  most  important  of  these  programs  are: 

1.  Program  S  provides  the  pressure  amplitude  and  phase  at  one  specific  point  in 
the  field.  It  must  be  used  with  the  "Record  on"  command  of  the  IBM  30-33  for 
sending  the  results  in  a  separate  file. 

2.  Program  SDF  is  a  version  of  the  DEME  program  which  takes  the  input  data 
from  an  input  file  and  transfers  the  output  to  another  file.  With  this  technique  the 
EASYPLOT  program  can  be  used  to  produce  multiple  graphs. 

3.  Program  PREC  uses  double  precision  to  reduce  round-off  errors.  The  outputs 
of  PREC  can  be  plotted  with  any  other  graphics  software  except  DISSPLA.  The 
DISSPLA  package  has  some  difficulty  when  applied  to  the  output  of  double 
precision  fortran  programs. 

4.  Program  SLAl  generates  graphs  without  the  slope  normalization. 

B.  PROGRAM  3LS4 

This  program  gives  alternate  forms  of  data  presentation.  It  uses  the  IBM  3033 
and  the  DISSPLA  package  to  generate  a  3-D  representation  of  the  pressure 
amplitude  distribution  on  the  bottom  or  along  any  plane  (defined  by  the  D)  within 
the  wedge.  It  also  provides  a  contour  presentation  of  the  same  data. 

Because  of  the  large  volume  of  numerical  output,  3LS4  is  ineffective  fo: 
recognizing  details. 

An  important  limitation  of  the  program  is  that  the  array  containing  the  data 
must  be  initially  dimensionalised  to  the  exact  size  of  the  data  file  array.  This 
information  is  normally  passed  from  the  data  file  during  the  program. 
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IV.  DISCUSSION  OF  THE  RESULTS 


The  DEME  and  3LS4  programs  have  been  created  using  the  same  model. 
Although  they  provide  different  output  graphs  they  must  agree  at  every  point.  The 
3LS4  provides  the  pressure  amplitude  and  phase  on  a  plane  defined  by  the  receiver 
angle  (D).  The  DEME  program  provides  the  pressure  distribution  and  phase  along 
a  vertical  line  anywhere  within  the  wedge.  The  programs  were  checked  for  both 
slow  and  fast  bottoms.  For  the  same  input  data,  they  gave  the  same  outputs,  as 
showm  in  Tables  1  and  2. 

For  a  slow  bottom  (ci  >  ci),  three  different  types  of  curves  of  pressure  amplitude  as 
a  function  of  depth  have  been  observed.  The  same  types  of  the  curves  have  been 
noticed  by  Kaswandi  [Ref  17]  and  are  shown  in  Figure  6. 

When  the  source  distance  (Rl)  is  much  smaller  than  the  characteristic  distance 
(Xj.  the  sound  pressure  is  linearly  dependent  on  receiver  angle  with  the  surface 
pressure  equal  to  zero  and  the  bottom  pressure  a  maximum  (Rl  =  0.1)  When  Rl 
approximates  X  ,  the  curve  develops  considerable  curvature  (Rl  =  0.7). 

When  the  source  distance  (Rl)  is  slightly  different  from  the  characteristic 
distance  X  (Rl  =  1.1),  then  a  maximum  and  a  minimum  appear  between  the 
surface  and  the  bottom. 

When  the  source  distance  (Rl)  is  greater  than  the  characteristic  distance  (X) 
(Rl  =  4),  then  the  sound  pressure  has  a  maximum  between  the  surface  and  the 
bottom.  No  minimum  appears  in  this  kind  of  the  curves.  When  the  value  of  Rl  is 
much  greater  than  X  ,  several  minima  and  maxima  appear  in  the  output  curve 
(Rl=  12). 
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Figure  7  shows  the  development  of  the  second  type  of  curve  (R1  =  1.0)  a?  the 
receiver  distance  is  varied. 

For  the  same  input  data  (  B  =  10,  G  =  5,  =  0-9),  the  fast  bottom  gave 

entirely  different  curves  as  shown  in  Figure  8.  When  the  source  distance  (Rl)  is 
less  than  three  times  the  characteristic  distance  (X),  only  the  first  mode  is  received 
everywhere  within  the  wedge  (Rl  =  0.3,  1.1,  1.3).  When  Rl  >  3,  higher  modes  are 
excited  and  the  curves  are  the  superposition  of  several  modes  (Rl  =  4). 

A.  TRANSITION  POINT 

For  the  slow  bottom,  when  the  source  distance  (Rl)  is  about  the  same  as  the 
characteristic  distance  (X),  the  curve  has  a  maximum  and  a  minimum.  The 
minimum  is  close  to  the  bottom.  A  “transition  point"  has  been  defined  by 
Kaswandi  [Ref  17]  as  the  receiver  distance,  for  which  the  pressure  above  a  minimum 
extrapolates  to  zero  on  the  bottom. 

Three  examples  are  shown  in  the  Figure  9  (R2  =  4.6,  6.4,  32).  In  all  of  them 
the  pressure  above  the  minimum  can  be  extrapolated  to  zero  pressure  point  at  the 
bottom.  The  transition  point  phenomenon  happens  for  source  distances  1<R1<1.5  . 

For  1  <  Rl  <  1.5  the  shape  of  the  pressure  curves  as  the  receiver  distance 
(R2)  is  varied  can  be  explained  as  follows:  as  the  receiver  is  moved  from  (R2  =  1.5 1 
to  (R2  =  30)  the  minimum  pressure  amplitude  decreases  in  magnitude  and  the 
first  transition  point  is  observed  (Fig  7).  With  further  increasing  of  R2.  the 
minimum  pressure  amplitude  of  the  curves  decreases  to  a  minimum  value,  and  then 
increases  arriving  at  a  second  transition  point.  For  the  receiver  distances  bet..een 
the  first  and  the  second  transition  point,  the  extrapolation  to  zero  pressure  was 
observed  to  be  above  the  bottom  (Figure  7).  With  further  increasing  of  R2,  a  third 
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transition  point  tries  to  form  but,  since  the  R2  is  so  large  round-off  error  appears 
in  the  program  output  curve  and  the  curve  characteristics  can  not  be  easily 
recognized.  In  the  figure  10,  the  locus  of  the  two  delectable  transition  points  are 
shown  as  a  function  of  R2  and  the  shore  distance  (YO).  Although  it  is  not  indicated 
in  the  figure  10  the  points  lie  on  circular  arcs  centered  at  the  source.  This  may 
indicates  that  the  transition  point  obeys  the  cylindrical  spreading  law. 

The  locus  of  the  first  and  second  transition  points  form  curves  which  can  be 
approximated  by  arcs  of  circles  centered  at  the  source  (figure  10) 

The  transition  point  as  a  function  of  the  source  angle  (G)  was  studied  using  the 
output  data  shown  in  Table  3.  As  the  source  angle  changes,  the  transition  point 
moves  to  a  different  receiver  distance.  Two  different  wedge  angles  were  observed. 
The  greater  the  source  angle  (G)  ,  the  smaller  the  receiver  distance  where  the 
transition  point  occurs. 

A  possible  explanation  is  that,  by  increasing  the  G  the  source  is  placed  far  from 
the  angle  value  (D),  where  the  minimum  happens,  so  the  absolute  pressure  values 
become  smaller  as  the  source  angle  increases.  In  figures  10a  and  10b  the  pressure 
distribution  of  the  bottom  is  shown  for  G  =  2.5®.  In  figures  10c  and  lOd  the  same 
case  is  represented  with  G  =  7.5®.  These  figures  are  generated  with  the  program 
3LS4  [Ref  18]. 

In  the  fast  bottom  case  no  transition  point  was  observed.  In  Figure  11  the  fast 
bottom  case  is  shown  when  the  shore  distance  (YO)  varies.  This  is  the  case  of 
downslope  propagation  in  a  wedge-shaped  fluid  overlying  a  fast  bottom.  At  the 
receiver's  position  (R2  =  32),  as  the  shore  distance  is  increasing  (  0  <  YO  <  80). 
only  a  curve  like  the  first  mode  appears  and  no  minimum  is  observed. 


16 


B.  MORE  CASES 

For  the  fast  bottom,  when  the  source  distance  (Rl)  is  bigger  than  3,  the  higher 
modes  start  to  appear.  In  Figure  12,  Rl  =  6  .  As  explained  by  Kaswandi  [Ref  17] 
when  Rl  =  6  the  first  three  modes  are  present.  By  increasing  the  source  angle  (G) 
the  output  curve  changes  because  at  some  source  positions  some  modes  are  reduced 
and  some  are  amplified. 

Some  problems  appear  when  the  source  is  very  close  to  the  apex  (Rl  =0.1),  as 
in  Figure  13.  The  program  seems  to  fail  close  to  the  bottom.  It  may  be  a  round-off 
error  of  the  computer  because  the  absolute  pressure  values  are  very  small.  The 
round-off  error  is  clearly  shown  in  Figure  14,  where  it  starts  to  appear  when  the 
source-receiver  distance  becomes  too  far.  An  attempt  to  eliminate  this  problem 
was  made  in  the  PREC  program  which  was  a  double  precision  version  of  the  DEME 
program.  Also  another  try  was  made  to  run  the  DEME  program  in  a  386  IBM 
compatible  personal  computer  with  16-bit  memory.  Both  of  the  tries  gave 
smoother  curves  but  they  did  not  solve  the  problem. 

C.  SHORE  DISTANCE  AS  A  VARIABLE 

Until  now  most  studies  were  for  on  axis  receiver  positions  (YO  =  0).  In  the 
following  discussion  the  shore  distance  (YO)  is  the  variable.  In  slow  bottom  three 
different  source  distances  (Rl)  were  studied  (Rl  =  0.5  ,  Rl  =  1  and  Rl  =  5). 

The  first  case  is  shown  in  Figure  15,  where  the  source  is  close  to  the  apex  (Rl= 
0.5).  At  receiver  position  R2  =  4  ,  only  the  linear  pressure  curve  is  observed  at  all 
values  of  YO  (0  <  YO  <  60). The  effect  of  larger  values  of  source  and  receiver 
distance  (Rl=  1,  R2=  6)  is  shown  in  Figure  16.  As  YO  increases,  the  curve  forms  a 
transition  point  at  YO  =  3  and  then  becomes  a  straight  line  at  YO  =  10. 
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Figure  17  illustrates  the  case  when  the  source  distance  is  greater  than  the 
characteristic  distance  (  R1  =  5).  Curves  for  four  different  values  of  YO  are  shown. 

The  same  cases  were  studied  also  for  the  fast  bottom.  The  outputs  are  shown  in 
Figure  18  when  the  source  distance  is  close  to  the  apex(Rl  =  0.5),  in  Figure  19 
where  R1  =  1  ,  and  in  Figure  20  where  the  R1  =  5.  What  is  clear  to  be  observed  in 
these  cases  is  the  uncertainties,  which  are  present  in  great  distances.  Also  when 
R1  >  3  the  higher  modes  start  to  appear. 

D.  ANGULAR  OUTPUTS 

Examination  of  the  contour  plots  revealed  a  general  hyperbolic  or  elliptical 
pattern,  therefore  another  try  was  made  to  understand  the  physical  structure  of  the 
acoustic  field  by  searching  at  some  particular  angle  off  the  axis  (  arctan[^’0/R2]  ). 
As  shown  in  Figure  21  output  data  were  calculated  every  10®  off  the  axis.  The 
acoustic  field  was  studied  for  two  different  source  positions,  for  a  small  source 
distance(Rl=0.5)  and  for  a  source  distance  greater  than  the  characteristic  distance 
(Rl=10). 

The  bottom  pressure  distribution  is  shown  in  Figures  22  through  25.  The  slow 
bottom  as  well  as  the  fast  bottom  has  been  studied.  From  the  contour  plots 
(Figures  23  and  25)  it  would  appear  that  large  values  of  the  axis  angle  were  more 
interesting,  because  the  isopressure  curves  seem  to  be  asymptotic  to  the  axis  angles. 

In  the  Figures  26  through  28  ,  the  pressure  distribution  along  three  different 
angles  are  shown  for  the  slow  bottom.  The  curves  are  unnormalized. 

The  fast  bottom  graphs  for  the  same  inputs  are  shown  in  Figures  29  through  34. 
Many  different  axis  were  studied  from  0®  to  90®.  All  the  curves  from  0®  up  to  70® 
seem  to  be  the  superposition  of  the  first  five  modes  (Figure  29,  30, and  31).  From 
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75®  axis  angle  and  for  big  receiver  distances  (R2)  the  graphs  begin  to  have  only  tlie 
first  mode.  That  means  that  the  higher  modes  have  been  cut  off.  Ai  80®  all  the 
outputs  appear  to  be  in  the  area  where  only  the  first  mode  is  received.  The  R2  in 
all  the  graphs  was  varied  between  10  and  60. 

The  slow  bottom  curves  for  any  angle  off  the  axis  were  observed  to  have 
uniform  vertical  pressure  profiles  as  showm  in  Figures  26  through  28.  The  fast 
bottom  curves  for  angles  up  to  80®  were  observed  to  have  complicated  profiles  as 
shown  in  Figures  29  through  32. 

D.  THE  OTHER  PARAMETERS 

Variations  of  the  other  parameters  such  as  the  velocity  ratio  Ci/c2  and  the 
density  ratio  Pj/P2  studied  by  Kaswandi  [Ref  17]  for  on  axis  measurements. 
The  results  showed  no  strong  variation  of  parameters  studied. 


19 


V.  CONCLUSIONS  AND  RECOMMENDATIONS 


In  this  research,  the  3-D  propagation  of  sound  has  been  studied  using  the  image 
theory.  Two  computer  programs  were  used  to  provide  the  pressure  amplitude 
distribution  and  phase  everywhere  within  a  wedge-shaped  fluid  overlying  a  fast  or  a 
slow  bottom. 

The  DEME  and  the  3LS4  programs  were  checked  for  the  same  input  data  and 
they  provided  the  same  outputs  for  the  same  acoustic  field  points.  After  this  check, 
it  was  concluded  that  the  programs  were  consistent  with  each  other. 

The  transition  point,  present  only  over  the  slow  bottom,  was  studied  at  various 
places  inside  the  wedge.  It  was  found  that  all  the  curves  could  be  approximated  by 
concentric  circular  arcs,  having  the  source  at  the  center. 

As  the  source  angle  (G)  increases,  it  was  concluded  that  the  transition  point 
distance  decreases  as  it  is  shown  in  Table  3  . 

When  the  source  position  is  1.1<  R1  <  1.5  then  two  transition  point  curves 
were  formed.  Sometimes  it  was  possible  to  have  also  a  third  one,  if  the  computer 
could  work  with  very  small  numbers  (16  or  32  bit  Bus). 

The  off  axis  data  were  taken  along  many  angular  axes.  The  slow  and  the  fast 
bottom  were  studied  for  the  same  angles  of  axis.  The  slow  bottom  data  showed  very 
small  variation  of  the  pressure  amplitude  values  as  the  receiver  distance  was 
increased.  The  fast  bottom  curves  showed  very  clearly  the  w'ave  guide  cutoff 
phenomenon.  In  ail  the  axis  angles  up  to  80®  the  curves  were  observed  to  be 
irregular  with  many  variations  and  unique  for  each  R2.  At  the  angle  of  80®  or 
more,  no  major  variations  were  observed  as  the  R2  was  increased.  This  observation 
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was  found  to  be  consistent  with  the  contour  plot  output  data  taken  from  the  3LS 1 
program. 

This  research  was  another  phenomenological  theoretical  step  toward  elucidating 
features  of  the  wedge  problem.  Some  experimental  data  must  be  compared  with  the 
two  programs  DEME  and  3LS4  to  have  some  confirmation  of  the  programs 
capabilities.  The  most  interesting  case  seemed  to  be  the  slow  bottom  case. 

Current  thoughts  [Ref  23]  suggest  that  these  features  can  be  studied  from 
another  point  of  view  based  on  a  collection  of  dipoles.  The  source  and  its  first 
image  form  a  dipole  as  do  other  associated  pairs  of  images  (Figure  1).  If  all  these 
pairs  of  images  are  considered,  then  the  field  can  be  considered  to  be  formed  from 
the  near  cancellation  of  equivalent  dipole  pairs  lying  at  slightly  different  distances 
from  the  receiver. 
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APPENDIX  A 


TABLE  1 

CONSISTENCY  OF  THE  DEME  AND  3L.S4  PROGRAMS  (SLOW  BOTTOM^ 
B  =  10®,  G  =  5®,  YO  =  50,  R1  =  40,  R2  =  10,  /)l/p2  =  0.9,  cl/c2  =  2 


RECEIVER 

ANGLE  D 

OUTCOMING 

DEME 

PRESSURE  AMPLITUDE 

3LS4 

00 

0.5073 

0.50728 

01 

0.8683 

0.86S28 

02 

1.4251 

1.42509 

03 

1.8969 

1.89691 

07 

1.8491 

1.84912 

10 

0.0000 

0.00000 
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TABLE  2 

CONSISTENCY  OF  THE  DEME  AND  3LS4  PROGRAMS  (FAST  BOTTON! ' 
B  =  10°,  G  =  5°,  R1  =  40,  R2  =  35,  =  0.9,  =  0.9 


REC.ANGLE  SHORE  DIST.  PRESSURE  AMPLITUDE  OUTPUT 


D 

YO 

DEME 

3LS4 

00 

00 

9.5446 

9.54476 

00 

30 

2.5776 

2.57740 

00 

50 

3.2210 

3.22089 
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TABLE  3 


TRANSITION  POINT  AS  A  FUNCTION  OF  SOURCE  ANGLE  (G) 

P^IP2  ~  ~ 


WEDGE 

ANGLE  (B) 

SOURCE 

ANGLE (G) 

SOURCE 

DIST.(Rl) 

RECEIVER 

DIST.(R2) 

10° 

2.5° 

1.0 

33.0 

10° 

5.0° 

1.0 

24.0 

10° 

7.5° 

1.0 

22.0 

o 

o 

1.5° 

1.3 

5.9 

o 

o 

3.0° 

1.3 

5.7 

06° 

4.5° 

1.3 

5.4 
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APPENDIX  B 


FIGURES 


Image  Structure  for  a  wedge-shaped  duct 
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Figure  3.  Rainge  and  Angle  at  Bottom  of  the  Image. 
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im 
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wn 


Figure  4.  Geometry  of  the  N***  Image  and  the  Receiver  projected  onto  a 

Plane 


28 


(03a)  310,VV  H3AI303a 


Pressure  Amplitude  at  R2s20  and  Y0=0  for  various  Rl  over  a 
Slow  Bottom.  P\Jp^  =  .9,  clarity  only  every  5*** 

computation  shown) 


KEC.ANGLE  VS.  NORM.  PRESSURE 


(33a)  310NV  a3AI333a 


Figure  7.  Pressure  Amplitude  at  Rlsl  for  Various  B2  over  a  Slow  Bottom 

=  ,9,  1.1  .(For  clarity  only  every  5^^  computation 

shown) 
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NORM.  RRKSSURE  AMPLITUDE 


NOKM.  VUKSSVKE  ANirLITUDE 


ft 


(oaa)  aiOiSiv  aaAiaoaa 


Figtire  9.  EIxamples  of  Transition  Points. 
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001 


Figure  10. 


08  09  O’V  02 

tea)  aoMvisia  aaAiaaaa 


Transition  Point  Curves 
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SLOW  BOTTOM 


^EC.  ANGLE  VS.  NORM.PRESSURE 


(oaa)  3i3is:v  aaAiaoaa 


Figure  11.  Preaeure  Amplitude  at  Rlsl.l  and  R2=32  for  Various  YO  over  a 
Fast  Bottom.  (For  clarity  only  every  computation  shown) 
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NORM.  PRESSURE  AMPLITUDE 


REC.ANGLE  VS.  PRESSURE 


I  •  •  »  . 

ST  cl  6  B  E  0 

(oaa)  310NV  a3Ai303a 


NORM.  PRESSURE  AMPLITUDE 


REC, ANGLE  VS.  NORM. PRESSURE 


NOW1.  PRESSURE  RHPLITUOE 


REC.ANGLE  VS.  NORM.  PRESSURE 


(oaa)  310MV  a3.\i333a 


Figure  14.  Prenure  Amplitude  at  RlsO.5,  Y0=0  for  Various  R2  over  a  Slow 
Bottom.  (For  clarity  only  every  5^**  computation  shown) 
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NORM.PRESSURE  AMPLITUDE 


REC.  ANGLE  VS.  NORM.  PRESSURE 


(03a)  330MV  a3AI333a 


Figure  15.  Pressure  Amplitude  at  RlsO.5,  R2=:4  for  Various  YO  Over  a 
Slow  Bottom.  (For  clarity  only  every  5^**  computation  shown) 
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NORM.PRESURE  AMPLITUDE 


REC.ANGLE  VS.  NORM.PRESSURE 


(03a)  aiOKV  H3AI333H 


Figtue  16.  Pressure  Amplitude  at  Rl=l,  R2s6  for  Various  YO  over  a  Slow 
Bottom.  (For  clarity  only  every  5^**  computation  shown) 
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NOKM.PRESSURE  AMPLITUDE 


(oaa)  313KV  a3Ai3D3a 


Figure  17.  Pressure  Amplitude  at  Rl=5,  R2=20  for  Various  YO  Over  a  Slow 
Bottom.  (For  clarity  only  every  computation  shown) 
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REC.ANGLE  VS.  PRESSURE 


(03a)  310NV  a3AI303a 


Figure  18.  Pressure  Amplitude  at  Rl=0.5,  R2=4  for  Various  YO  Over  a  Fast 
Bottom.  (For  clarity  only  every  5^**  computation  shown) 
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NORM.  PRESSURE  AMPLITUDE 


REC.ANGLE  VS  PRESSURE 


(93a)  339MV  a3AI39aa 


Figure  19. 


Pressure  Amplitude  at  Rl=l,  R2s:6  For  Various  YO  Over  a 
Fast  Bottom.  (For  clarity  only  every  5***  computation  shown) 


NORM.PRESSURE  AMPLITUDE 


REC.ANCLE  VS.  PRESSURE 


(030)  330NV  aaAi3Daa 


Figure  20.  Pressure  Amplitude  at  Rl=5,  R2=20  for  Various  YO  Over  a  Fast 
Bottom.  (For  clarity  only  every  5^^  computation  shown) 
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NOKM.PUESSURE  AMPLITUDE 


BOTTOM 


REC.  ANGLE  VS.  PRESSURE 


(03a)  330NV  asAoosa 


Figure  26.  Pressure  Amplitude  at  60°  off  Axis  (Rl=10)  for  Various  R2  over  a 
Slow  Bottom.  (For  clarity  only  every  5^**  computation  shown) 


5' 


0.3  0.6  0.9 

PRESSURK  AMPLITUDE 


0.0  0.3  0.0  0.9 

TRESSURE  AMPLITUDE 


REC.  ANGLE  VS.  PRESSURE 


c 


(03a)  330i\:v  a3.\i303a 


Figure  28.  Pressure  Amplitude  at  80^  off  Axis  for  Various  R2  (RlsslO)  over  a 
Slow  Bottom.  (For  clarity  only  every  5***  computation  shown) 
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0.0  0.2  0.4  O.G  0.0 

PRESSURE  AMPLITUDE 


REC.  ANGLE  VS.  PRESSURE 


(030)  300NV  a3AI303a 


Figure  30. 


Pressure  Amplitude  30^  of  Axis  for  Various  R2  (Rl=10)  over  a 
Fast  Bottom.  (For  clarity  only  every  5^**  computation  shown) 


2.0  3.0 

PRESSURE  AMPLITUDE 


REC.  ANGLE  VS,  PRESSURE 


OT  9  9  t  3  0 


(oaa)  aioNV  aaAiaoan 


Figure  31.  Pressure  Amplitude  at  60^  of  Axis  for  Various  R2  (Rl=10)  Over  a 
Fast  Bottom.  (For  clarity  only  every  5^  computation  shown) 
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0.0  1.0  2.0  3.0  4.0 

PRESSURE  AMPLITUDE 


REC.  ANGLE  VS.  PRESSURE 


(oaa)  310NV  aaAiaoaa 


Figure  32.  Prenure  Amplitude  at  75^  of  Axis  for  Various  R2  (Rl=10)  Over  a 
Fast  Bottom.  (For  clarity  only  every  5^^  computation  shown) 
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0.0  1.0  2.0  3.0  4.0 

PRESSURE  AMPLITUDE 


REC.  ANGLE  VS.  PRESSURE 


(oaa)  aioNv  aaAiaoaa 


Figure  33.  Pressure  Amplitude  at  77.5^  of  Axis  for  Various  R2  (RlslO)  Over 
a  Fast  Bottom.  (For  clarity  only  every  5^**  computation  shown) 
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00  I.o  2.0  3.0 

PRESSURE  AMPLITUDE 


REC.  ANGLE  VS.  PRESSURE 


(03a)  330K;V  H3AI303H 


Figure  34.  Pressure  Amplitude  at  80^  of  Axis  for  Various  R2  (Rl=10)  Over  a 
Fast  bottom.  (For  clarity  only  every  5'**  computation  shown) 
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0.0  0.3  0.0  0.9 

PRESSURE  AMPLITUDE 


APPENDIX  C 
DEME  PROGRAM 


This  program  calculates  the  pressure  amplitude  and  phase  everywhere  within 
a  wedge  shape  fluid  overlying  a  penetrable  bottom.  It  is  a  modified  version  of  the 
DSLOW  proCTam,  which*  was  written  by  CASWANDI.  It  has  been  modified  by 
LT.D.PALIATSOS  on  March  89.  This  program  is  working  for  fast  or  slow  bottoms. 
The  output  data  are  given  bv  a  table  or  by  a  graph  using  DISSPLA. 

»*7if**  'if  »i***^^**^^t|^^i:^*m** 

INTEGER  A,I,I1,M,N,S1,S2,N1,J,K 

REAL*4  B,CC,C2,D,D1,D2,G,PI,P1,P2,Q1,R1,R2,T, 

*  T4,T6,VV0,W1,Y0,Y1,Y2,Z1,Z2,Z3,Z4,Z5,Z6, 

*  T1{80),R8(80),R9(80),S(30),C(30),E(30), 

*  F(30),Y,Z,R3,AL,PZ(30),DZ(30),V,DX,XP,DD,PN(30) 

REAL*4  TQQ,TQQ1,TQQ2,TQQ3 

PI  =ACOS(-1.0D00) 

^  Hi  DC « « Ik  :tc  *  «  « iti « :tc «  « :(i  lit  itiili  m  « i|c  ^  Ik  i|l »  A  »  «  *  1)1  *  1*1 iti  :*[  !([  4c  k  *  «  «  *  *  itt  *  Ik  *  DC  Ik  X  X  X  *  K 

C  INPUT  PARAMETERS 

Qikik4cicxikx4cikikikikikikHcikik4c4cikHcxxxxikikHc4<ikik4cxikik4c4c4cx4cik4c4cikik4cikikikikikikikHcikikikxx4cx*ikik>kikxxxxx 

C  B  =  WEDGE  ANGLE  (DEG) 

C  G  =  SOURCE  ANGLE  (DEG) 

C  D  =  RECEIVER  ANGLE  (DEG) 

C  Nl=  #  OF  IMAGE  POINTS 
C  R1  =  SOURCE  DISTANCE  (IN  DUMP  DISTANCES) 

C  R2  =  RECEIVER  DISTANCE  (IN  DUMP  DISTANCES) 

C  YO  =  APEX  DISTANCE  (IN  DUMP  DISTANCES) 

C  D1  =  RHO  1/RHO  2 

C  CC  =  C-l/C-2 

C  AL  =  ALPHA/K2 

C  A  =  #  OF  RECEIVER  POSITIONS 

^ikikxikxXikxikHcikikDcDcDcDcxDcxHcxDcikDcXHcxxikxxxscikikikikDcikikHcikDcikHcikDiikDcikikikikikikikikikikKikikDcHcxxxKikiKx 

C  INITIAL  INPUT  RELATIONS 

^  X  X  X  « A  »  «  « 3K  3K  Jk  «  X  X  *  ^  «  4c  *  V  *  ^  lit  *  ](c  Jtc  X  *  Sic  Mc  *  X  X  X  X 

B  =  5.5 
G  =  2.75 
D  =  00.0 
R1  =  1.3 
R2  =  9.0 
YO  =  2.0 
D1  =  0.80 

C  DO  90  P  =  1,4 

C _ 

c  YOU  must  Give  the  Information  aBOut  T'He  BOTTom;  That 

C  MEANS  YOU  MUST  CHOOSE  SLOW  OR  FAST  BOTTOM  BY  GIVING 
C  SOME  VALUE  TO  THE  SPEED  RATIO  CC  =  C1/C2  .  ] 

C _ 

C  CC=C1/C2 
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CC  =1.1 
AL  =  0.0001 
A  =  10 

0  «  »  *  *  at  Hi  Di  H  *  m  Ik  iti  *  :ti  4c  :«l  Jk  *  «  »  * *  *  1)1 «  •  :ti «  DCDC  DI »  *  4i  He  *  m  Ik  » :|c  itt  *  DE  Ik  «  Hi  :«E  *  V  K  X  K  ai  ac »  ac  X  x  DC  X  x  X  X 

C  MAIN  PROGRAM 

Q  X  Hi  Ik  Hi  Ik  *  Hi  HI  Hi  Hi  Hi  Ik  Ik  Hi  Hi  He  Ik  HI  HI  *  Ik  Ik  Ik  H<  *  *  Hi  Hi  K  *  Hi  «  Hi  X  Hi  Hi  Ik  Hi  Ik  )k  iti  HI  Ik  Hi  <1 W  Ik  He  Hi  H;  Hi  Ik  HI  He  He  H;  He  H;  *  Ik  Hi  X  X  X  X  X  X  X  X 

K  =  0 

N1  =  INT(180./B) 

T6  =  180./PI 
B  =  B/T6 
G  =  G/T6 
C2  =  00**2 

D2  =  (Y0*Y0)+(R1*R1)+(R2*R2) 

R3  =  2  *Rl*R2 

0  T4  =  PI/(2*TAN(A0OS(l/00))*TAN(B)) 

TQQ  =  TAN(B) 

C - 

0  DEOISION  ABOUT  SLOW  OR  FAST  BOTTOM 

C - 

IF  (OO.LT.l)  THEN 
TQQl  =A0OS(00) 

TQQ2  =  SIN(TQQl) 

0  T4  =  PI/(2*SIN(AC0S(00))*TAN(B))  *****  FAST  OASE 
ELSE 

TQQl  =  A00S(l/00) 

TQQ2  =  TAN(TQQl) 

0  T4  =  PI/2*TAN(A0OS(l/00))*TAN(B)  ****  SLOW  OASE 
ENDIF 

0 - 

TQQ3  =  2.*TQQ2*TQQ 
T4  =  PI/TQQ3 

0  K1X=  PI/(2*TQQ*TQQ2) 

D2  =  Y0*Y0+RrRl+R2*R2 
Q1  =  1/DSQRT(2.0D00) 

0300  FORMATCNOP  ',7X,'THETA(N) 'JX.TMGE  SR  R8  ',7X, 

0  -  TMGE  SR  R9  ') 

0  WRITE  (6,300) 

800  FORMAT(5X,'REO.POS',5X,'  REO.ANGLE  ',5X,'PRES.AMPLITUDE  ',5X. 
*  ’PHASE  ANGLE  ',5X,'NORM.PRESS') 

WRITE  (6,800) 

DO  10  M  =  0,A 
D  =  M*B/A 
DD=B/A 
V  =2*6/10 

IF(D.LT.V)  GOTO  110 
IF(D.GE.V)  GOTO  120 
110  DX  =  D 

DO  15  J  =  1,10 
D  =  DX+(J-1)*B/(10*A) 

DD  =  B/(10*A) 

120  SI  =  1.0 
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DO  20  N  =  1,N1 
IF(Sl.GT.O)  Tl(N)=(N-l)*B+G 
IF{S1.LT.0)  T1(N)=N*B-G 
SI  =  -  SI 

RS(N)  =  SQRT(D2-R3*C0S(T1(N)-D)) 
R9(N)  =  SQRT(D2-R3*C0S(T1(N)+D)) 
C310  FORMAT(3X,I2  ,5X,F5.2,4X,F6.4,6X,F6.4) 
C  WRITE  (6,310)  N,T1(N),R8(N),R9(N) 

20  CONTINUE 


PI  =  0.0 
P2  =  0.0 
DO  30  N  =  1,N1 
S2  =  (-l)**(INT(N/2)) 

W1  =  2=^C2*AL 
II  =  INT((N-l)/2) 

DO  40  I  =  1,11 

S(I)  =  ABS(R1*SIN(T1(N)-2*I*B) 

*  +R2*SIN(2*rB-D))/R8(N) 

C(I)  =  SQRT(l-(S(I)*Sa))) 

T  =  S(I)/D1 

WO  =  (-C2+(C(I)*C(I))) 

Y  =  SQRT((W0*W0)+(W1*W1)) 

Z  =ABS(W0) 

IF(Y.LE.Z)  Y  =  Z 
Y1  =  Q1*SQRT(Y+W0) 

Y2  =  -<J1*SQRT(Y-W0) 

Z1 =  T-Y2 
Z2  =  -Y1 

C  Z12=  CMPLX(Z1R.Z2I) 

Z3  =  Z1/(Z1’‘Z1+Z2*Z2) 

Z4  =-Z2/(Zl*Zl+Z2*Z2) 

Z1  =  T+Y2 
Z2  =  Y1 

C  Z13=  CMPLX(Z3R.Z4I) 

Z5  =  ZrZ3-Z2’'Z4 
Z6  =  Z1*Z4+Z2*Z3 
C  Z14=  CMPLX(Z5R,Z6I) 

E(I)  =  Z5 
F(I)  =  Z6 

C400  FORMATC  NOP  ',5X,'I  ',5X,'ARSIN(S(I)*T6)',5X,'E(I)', 

C  *  5X,'F(I)',5X,  '  '  '  '  / 

C  *  '  bSQRT(E(I)’‘E(I)+F(I)*F(I))',5X,'T6*ATAN(F(I)/E(I))’) 

C  WRITE  (6,400) 

C420  FORMAT(3X,I2.6X,I2,5X,F6.4,10X,F6.4,3X,F6.4,3X,F6.4,10X.F7.4) 
C  WRITE  (6,420)  N,I,  ASIN(S(I))*T6,E(I),F(I), 

^  SQRT(E(I)*E(I)+F(I)-F(I)),ATAK(F(I)/E(I)) 

40  CONIIinUE 
Z1  =  0 
Z2  =  0 
Z3  =  0 
Z4  =  0 
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Z5=  1 
Z6  =  0 

IF(N.LE.2.00)  GOTO  50 
DO  45  I  =  1,  II 
Z1  =  E(I) 

Z2  =  F(I) 

Z3  =  Z5 
Z4  =  Z6 

Z5  =  Z1*Z3-Z2*Z4 
Z6  =  Z1*Z4+Z2*Z3 
45  CONTINUE 
50  Z1  =  Z5 
Z2  =  Z6 
T  =T4*R8(N) 

Z3  =  COS(T) 

Z4  =  -SIN(T) 

Z5  =  Z1*Z3-Z2*Z4 
Z6  =  ZrZ4+Z2*Z3 
PI  =  P1+S2*Z5/R8(N) 

P2  =  P2+S2*Z6/R8(N) 

C500  FORMATC  NO.OF  I.P  ',3X,'RE(REFL)=  ',3X,TM(REFL)=  ') 

C  WRITE  (6,500) 

C510  FORMAT(3X,I2,6X,F6.4,6X,F6.4) 

C  WRTTF  (f\  ^9*70 

C600  FORMATC  I',5X,'ARCSfN(S(I))*T6=',3X,'E(I)',5X,T(I)’,5X.'EF-, 
C  *  5X,'  AN2  ') 

C  WRITE  (6,600) 

11=11+1 
DO  60  I  =  1,  II 

S(I)  =  ABS(RrSIN(Tl(N)-2*(I-l)*B) 

*  +  R2*SIN(2*(I-iyB+D))/R9(N) 

C(I)  =  SQRT(1-S(I)*S(I)) 

T  =  S(I)/D1 

WO  =  -C2+C(I)*C(I) 

Y  =  SQRT((W0*W0)+(W1*W1)) 

Z  =  ABS(WO) 

IF(Y.LE.Z)  Y  =  Z 
Y1  =  QrSQRT(Y+W0) 

Y2  =  -QI*SQRT(Y-W0) 

Z1  =  T-Y2 
Z2  =  -Y1 

Z3  =  Zl/(ZrZl+Z2*Z2) 

Z4  =  -  Z2/(ZrZl+Z2*Z2) 

Z1  =  T+Y2 
Z2  =  Y1 

Z5  =  ZI*Z3-Z2’^Z4 
Z6  =  Z1*Z4+Z2*Z3 
E(I)  =  Z5 
F(I)  =  Z6 

AN1=  ASIN(S(I))*T6 

EF  =  SQRT(E(I)*E(I)+F(I)*F(I)) 
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AN2=  ATAN(F(I)/E(I))-T6 

C610  FORMAT(2X,I2,5X,F6.4,5X,F6.4,5X,F6.4,5X,F6.4,5X.F7.4) 

C  WRITE  (6,610)  I,AN1,E(I),F(I),EF,AN2 
60  CONTINUE 
Z1  =  0 
Z2  =  0 
Z3  =  0 
Z4  =  0 
Z5  =  1 
Z6  =  0 

DO  80  I  =  1,11 
Z1  =  E(I) 

Z2  =  F(I) 

Z3  =  Z5 
Z4  =  Z6 

Z5  =  Z1*Z3-Z2*Z4 
Z6  =  Z1*Z4+Z2*Z3 
80  CONTINUE 
Z1  =  Z5 
Z2  =  Z6 
T  =T4*R9(N) 

Z3  =  COS(T) 

Z4  =  -SIN(T) 

Z5  =  ZrZ3-Z2‘Z4 
Z6  =  Z1*Z44-Z2*Z3 
PI  =  P1+S2*Z5/R9(N) 

P2  =  P2+S2*Z6/R9(N) 

C700  FORMATC  LOWER  PATH  NO=',3X.'RE(REFL)=  ’.3X.'IM(REFL)=  ') 
C  WRITE  (  6,700) 

C710  F0RMAT(6X,I2,12X,F6.4,12X,F6.4) 

C  WRITE  (6,710)  N,S2*Z1,S2*Z2 
30  CONTINUE 
K=K-I-1 
DZ(K)=D*T6 

PZ(K)=SQRT(PrPl+P2*P2)*Rl 
C  WRITE  (6,810)  K,DZ(K),PZ(K),ATAN(P^P1),PN(K) 

C810  FORMAT(6X,I3,11X,F5.2,12X,F7.4,12X,F7.4,10X,F7.4) 

15  CONTINUE 
10  CONTINUE 
DO  31  L=l,29 
XP=PZ(28) 

PN(L)=PZ(L)/XP 

WRITE  (6,811)  L,DZ(L),PZ(L).ATAN(P2/P1),PN(L) 

81 1  FORMAT(6X,I3,l  1X,F5.2,12X,F7.4,12X,F7.4,10X,F7.4) 

31  CONTINUE 

C250  FORMATC  WEDGE  ANGLE  =  •,F5.2, IX,' SOURCE  ANGLE  ',F4. 

C  *  '  SOURCE  DISTANCE=',F4.2, IX, 'RECEIVER  DISTANCE=  ',F4. 

C  *  '  SHORE  DISTANCE=  ',F4.2) 
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(M  iM 


C  WRITE  (6,250)  B*T6,G*T6.R1,R2,Y0 

C270  FORMATC  RH01/RH02=',F5.2,5X,’  Cl/C2=‘.Fo.2.5X.' 

A1PHA/K2=',F8.4) 

C  WRITE  (6,270)  D1,CC,AL 
C  YO  =  YO+2 

C  R2  =  R2+2 

C  R1  =  Rl+2 

90  CONTINUE 
B  =  B*T6 
G  =  G*T6 
C  STOP 
C  END 

Q  Itciti  *  *  «  *  *  *  A  *  *  *  *  Hi  *  itl  *  IK  DC  *  «  *  « iti  *  *  *  4c «  «  *  *  «  *1  *  1)1  *  » 1|I  I|I  *  *  *  1*: «  *  Dt  He  *  iti  :tc  *  X  W  W  »  «  x  V  *  V  »  X  X  xx 

C  A  PROGRAM  FOR  PLOTTING  BY  TEK618  OR  SHERPA 

^HiX*HiHiHiHiHiHiHiHcXXX4cHiHiHcHiXHcXHiXHcHixHcXXXHiHcHiHcXHiHiHiHiHiHiHiHiHiXHi»XXXHixHiXHcXXXx*xxxXxXxx 

CALL  MEDBUF 
CALL  TEK618 

C  CALL  SHERPA('SLOWDEMl','A',3) 

CALL  NOBRDR 
CALL  PAGE(15.,12.) 

C  CALL  PAGE(8.5,11.) 

CALL  HWROT('AUTO') 

C  CALL  AREA2D(09.,5,) 

CALL  AREA2D(09.,7.) 

C  CALL  AREA2D(8.5,6) 

C  CALL  AREA2D(6.0,8.5) 

CALL  HEIGHT(.2) 

CALL  XNAME('NORM.  PRESSURE  AMPLITUDES'.24) 

CALL  YNAME('RECEIVER  ANGLE(DEG)S',19) 

CALL  YTICKS(5) 

CALL  XTICKS(5) 

CALL  GRAF(0.,2.0  ,10.0,0., 1.0, 6.0) 

C  CALL  GRAF(0.,1.0  ,05.0,0.,! .0,06.0) 

CALL  DOT 
C  CALL  GRID(2.2) 

CALL  HEADIN('REC. ANGLE  VS.  NORM.PRESSURES'.-l 00, 1.8.1) 

CALL  MESSAG('WEDGE  ANGLE  =$',100,8.,7.) 

CALL  REALNO(B,2,10.2,7.) 

CALL  MESSAG('RH01/RH02  =$',100,8.,6.5) 

CALL  REALNOfDl, 2, 10.2,6.5) 

CALL  MESSAG('C1/C2  =$',100,8., 6.0) 

CALL  REALNO(CC,2,10.2,6.) 

CALL  MESSAGf’SOURCE  ANGLE=  $',100,8.,5.5) 

CALL  REALNO(G, 2,10.2,5.5) 

CALL  MESSAGf'SOURCE  DIST.=  $',100,8.,5.) 

CALL  REALNO(R1,2,10.2,5.) 

CALL  MESSAGf'REC.DIST.  =  $',100,8.,4.5) 

CALL  REALNO(R2, 2, 10,2,4.5) 

CALL  MESSAGf'SHORE  DIST.=  $',100,8.,4.) 

CALL  REALNO(Y0,2,10.2,4.) 

CALL  MESSAG('PRESS.AMPL.(X=1)  =  S',100,7.0.3.5) 
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CALL  REALNO(XP, 105,10.2,3.5) 
CALL  RESET('ALL’) 

CALL  PARA3 
CALL  NOCHEK 
CALL  CURVE(PN,DZ,29,1) 
CALL  ENDPMO) 

CALL  DONEPi 

STOP 

END 
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